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abstract 

In this paper we shall introduce a lattice model of unconventional superconductors (SC) like d-wave 
SC in order to study quantum phase transition at vanishing temperature (T). Finite- T counterpart 
of the present model was proposed previously with which SC phase transition at finite T was in- 
vestigated. The present model is a noncompact U(l) lattice-gauge-Higgs model in which the Higgs 
boson, the Cooper-pair field, is put on lattice links in order to describe d-wave SC. We first derive 
the model from a microscopic Hamiltonian in the path-integral formalism and then study its phase 
structure by means of the Monte Carlo simulations. We calculate the specific heat, monopole densi- 
ties and the magnetic penetration depth (the gauge-boson mass) . We verified that the model exhibits 
a second-order phase transition from normal to SC phases. Behavior of the magnetic penetration 
depth is compared with that obtained in the previous analytical calculation using XY model in four 
dimensions. Besides the normal to SC phase transition, we also found that another second-order 
phase transition takes place within the SC phase in the present model. We discuss physical meaning 
of that phase transition. 
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1 Introduction 



Ginzburg-Landau (GL) theory plays a very important role in study of the superconducting (SC) 
phase transition. In particular GL models defined on a lattice have been extensively studied in order 
to investigate effects of topological excitations like vortices nonperturbatively. These models are 
nothing but gauge-Higgs models on a lattice that also play an important role in elementary particle 
physics. 

At present many materials that have unconventional SC phase have been discovered [1]. Among 
them, some of the materials, including the high-temperature (T) cuprates, the rare-earth heavy- 
fermion materials, etc, exhibit d-wave SC. Order parameter (the Cooper-pair wave function) of the 
d-wave SC has nodes. Therefore in order to formulate a GL theory of the d-wave SC, the Cooper-pair 
field, i.e., the Higgs field, must be put on links of the lattice instead of lattice sites. 

In the previous papers[2, 3], we introduced a new type of lattice gauge-Higgs models that describe 
d-wave SC phase transitions. We call these models U-V models where the [/-field refers to the 
ordinary gauge field whereas the V- field is the Higgs field defined on links. In Rcf. [2], a compact- 
U(l)-gauge version was studied from the viewpoint of the t-J model that is a canonical model of the 
high-T c cuprates. In the slave-particle representation, the [/-field is an emergent compact U(l) gauge 
boson and the V-field is a spinon-pair field for the resonating-valence-bond (RVB) order. By means 
of Monte Carlo (MC) simulations, we found that there appears a phase transition line as the coupling 
between the U and V-fields is increased. We measured instanton densities and expectation value of 
the Wilson loops, and verified that the observed phase transition is the one from a confinement to 
Higgs phases. There the confinement phase corresponds to electron phase, whereas the Higgs phase 
to fractionalizcd spin-gap phase. Numerical results show that the phase transition is of first order. 

In Ref. [3] , on the other hand, we studied finite-T phase transition of unconventional SC by using 
the U-V model. In this approach, the [/-field is a noncompact U(l) gauge field corresponding to 
the usual electromagnetic field, whereas V is the electron Cooper-pair field. There we considered 
the London limit of V, i.e., V x i SU(1). Gauge-invariant gauge-boson mass (the inverse correlation 
length of the magnetic field) is an order parameter for the normal to SC phase transition. We 
observed the phase transition by means of the MC simulations and found that the transition is of 
first order. However as the anisotropy of the layered structure is increased, signal of the first-order 
phase transition is getting weaker, e.g., the discontinuity of the gauge-boson mass at the transition 
point is getting smaller. 

In the present paper, we shall continue to study critical properties of the unconventinal SC by 
using the U-V model. In particular we focus on the quantum phase transition (QPT) at T = 0. 
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In Section 2, we shall derive the GL model for the QPT starting from a microscopic Hamiltonian 
of the Cooper pair. In Section 3, we shall show the numerical calculations of the U-V model in 
(2 + 1) dimensions ((2 + 1)D). We calculated the "internal energy", "specific heat", "monopole" 
density of the U and U-fields, and the magnetic penetration depth (the inverse gauge-boson mass). 
We found that the model exhibits a second-order phase transition from the normal to SC phases at 
certain critical value of the parameter g c , which controls quantum fluctuations. (Larger g suppresses 
quantum fluctuations more.) In Section 4, we report the result of the MC simulations for the U-V 
model in (3 + 1)D. The reason why we study both the (2 + 1) and (3 + 1)D models is that most of the 
materials of the d-wave SC are quasi-2D and have a layered structure. We found that the (3 + 1)D 
model also exhibits a second-order SC phase transition. However besides that transition, the model 
has another phase transition at larger g' c > g c within the SC phase. This phase transition is also of 
second-order. We shall discuss its physical meaning. Section 5 is devoted for discussion. 



2 Ginzburg-Landau Model on a lattice 

In this section, we shall derive the lattice GL model from the microscopic Hamiltonian by using the 
path-integral formalism. To this end, let us first introduce a Cooper-pair field on links (x, i) of 
d-dimensional hypercubic lattice with lattice spacing a^{~ 1), where x denotes the lattice site and i 
is the direction index i = 1, • • • , d. In terms of the electron annihilation operator il> XlCr , where a =f, | 
is the spin index, 

®x,i = i> x ,iipx+i,l ~ V'x.iV'x-KT- ( 2J ) 

In the d-wave SC's, amplitude of on-site Cooper pair is vanishingly small because of, e.g., the strong 
on-site Coulomb repulsion. Furthermore, for example, in the d x 2_ y 2-w&ve SC state, {<& X: i ■ <& x ,2) < 0, 
i.e., the Cooper-pair field changes its sign under a tt/2 rotation in the 1 — 2 plane of the real space. 
One should also notice that SC phase transition in the underdoped region of the high-T c materials 
can be understood as a Bose-Einstein condensation of a bosonic bound state of a pair of charge 
carriers. From these observations, we propose the following quantum Hamiltonian for the d-w&ve 
SC's, 

H = t {^x,iU* +i ^* x+j ,iKo + H ' c -) - C E ^:.,'K, J 'K.,.,'K, + H.c.) 

x,i=jtj x,i=£j 

+D (fxM^Ul+j.V:,! + H - c -) + E V x , y {N x ,i - Po )(N yj - po), (2.2) 

x,i=£j x,y,i,j 

where t, C, D are parameters, the electromagnetic (EM) gauge field U x ^ is related to the EM vector 
potential A cm as U x ^ — e 1 ^- 1 , A x ^ = J^ +l A em d£ and N X:i is the number operator of the Cooper- 
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pair field, N Xi i = $* ^ x ,i- Po denotes an uniform background positive charge and V x . y is a Coulomb 
potential including long-range interactions like, 

{e 2 Vn , x = v 

(2.3) 
e 2 /\x-y\ 1 |a:-y|»l 

for d = 3. As wc arc interested in the long-range interaction in Eq.(2.3) in the later discussion, 
we have ignored the dependence on direction of the Cooper pairs in the Coulomb potential. The 
i-term in Eq.(2.2) is noting but the hopping of the Cooper pair and the C-term is a plaquette 
term of the link field $> x j that stabilizes the phase fluctuations of & Xi i. On the other hand, the 
Z)-term determines relative phase of adjacent $'s. The Hamiltonian (2.2) might be derived from 
a canonical model of the high-T c cuprates like the t-J model by integrating out electrons just like 
Gorikov's derivation of the GL Hamiltonian for the conventional s-wave SC [4] [5]. It is obvious that 
the Hamiltonian H in Eq.(2.2) is invariant under the following local gauge transformation, 

U Xti -e^+'t^.ie-*-, $ x ,<-»e- <w "+'$ x , i e- ia '". (2.4) 

In the rest of this section, we shall derive the GL model from the Hamiltonian (2.2). If the reader 
is not interested in the derivation, please go to the result Eq.(2.27) directly. 

We employ the path-integral formalism by introducing imaginary time r. Then the partition 
function Z is given by 

Z = J [d<S>dA}e- s , (2.5) 

with the following action S 

S = f dT [Y,®U d ^^ + H ^iU)\ +S cm (A), (2.6) 

x,i 

where 

S cm {A) oc J drJ2F 2 (x) = J dTj2(A Xjl - A x+j4 + A x+iJ - A x .jf. 
Let us decompose $ x .i into radial and phase variables as follows, 

*x,i = (^ + ^x,i)e i0 - i . (2.7) 

By substituting Eq.(2.7) into Eq.(2.6), the first term of S becomes 

x,i 

= j dTY2i^ x ^^5 Px , l + 0{(8 Px ^) 2 ) 1 (2.8) 



4 



where we have used the periodic boundary condition on 5p x ^ and <f> Xt i (mod 2tt) in the r-direction, 
and <j> Xi i — d r (j) x .i. Similarly the terms in H Eq.(2.2) are expressed as 

H t ee -t Y^($ xA U: +iJ $l +jii U* J + c.c.) 

x,i=£j 

= -2tp ^ cos (0a:,i ~ <t>x+i,j ~ A xJ - A x+i j) + 0(5p x ^), 
x,i=jtj 

h c = c Y. : - ,|, :, ; • <••<••• 

x,i=£j 

= -2C>0 ^ COS (^,i ~ <j>x+i,j + 4>x+j.i - <f>x,j) + 0{8p x ,i), 
x,i=£j 

x,i^j 

= 2Dpl ^ C0S (^x,i - <Px+i,j - A x+j .i - A xJ ) + 0(Sp Xti ), 

x,i=£j 

Hn = Y V XiV (N Xii - p )(N y j - p ) 

x-,V,i,j 

= 1V .-//':"V',.,'V',.., • 0{{6p Xii f). (2.9) 

x,y,i,j 

From Eqs.(2.8) and (2.9), we can integrate over 8p x ^ in Eq.(2.5) and obtain, 

e' s * ee j[d5p]e-^-i dTHN 

= e -S dT T. x ,i,j K,l<f>x.ij>y,j _ (2.10) 

Furthermore, we can prove the following identity by straightforward calculation (see Appendix A), 

e- s * = J[dA ]e-f dTHA , (2.11) 

with 

H A = Jdk d 2 [e- 2 V(k) - l]- 1 i (-fc)io(fc) + ^ J2(4>x,i - eA xfi )(j> Xij - eA x . Q ), (2.12) 

x,i,j 

where V(k) is the Fourier transformation of V XlV and the scalar potential A x ^ is the field that 
mediates the Coulomb interaction (2.3). (Ao(k) is the Fourier transformation of A X fi.) 
Similarly we can show (see Appendix A), 

e -JdrH t = J[dV*dV}e--T dTHv , (2.13) 
\ E [j^e^V^ - V Xti e-^"* + c.c] , (2.14) 

x,ijtj 



where 
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with A = A x+ ij + A x j and J = 2tp . As V x+ j.i cx e~ %A+% ^ x -*, the field V Xi i can be regarded as the 
Cooper-pair field. 

Finally we shall integrate over <f) X:i . To this end, we consider the following integral 
e -fdrH' v = J ^ exp / J dT _ J_ J2(4> x , t - eA Xt0 )(4> xJ - eA xfi ) 

+(d - 1) ^(V^e"'*-' + c.c.)] ] • (2.15) 

x,i / 

From Eq.(2.15), it is obvious that requirement of the invariance of H' v under a "local gauge trans- 
formation" , 

4x,i^<t>x,i + ea x {r), A xfi -» A xfi + d x (r), V x>i -» e iea *^V Xti , (2.16) 
determines dependence on A x , in ffy. Then by simply putting A x $ — 0, we define 

(O) =Z^ 1 J[d(j>} e-pf^+'-'+'-'O, = J[d<j)} e -7*f dT ^ (2.17) 

It is not difficult to show, 1 

e 2 

{<f>x,i(n)4>v,ifo))o = Sxy lim £ ^ -e- e l ri - T2 l, (2.18) 
where e is the infrared cutoff. Then, 

( e m>.An)-4>.An))) = cxp [_ e 2| ri _ T2 |]. (2.19) 

From Eqs.(2.15) and (2.19), the lowest-order terms of the cumulant expansion are obtained as 
-i^Uo=o = (rf-l) 2 (/rfr 2 ^^(r)^(r 2 )e-^M-^(-))) o 

= (d - I) 2 1 rfT 2 \4, i (r)^(r 2 )e- e2 l— =1 

= A!]r|V^(T)| 2 -^|^(T)| 2 + ..., (2.20) 

x, i x,i 

where 

\ 1 = {d-lf J dT e -« a M=(d-l)2^, 

A 2 = (d - l) 2 y drr 2 e- e2 l T l = (d - l) 2 ^. (2.21) 

1 It should be remarked here that the Green function like {e l ^ a =. i ( T1 - > ~^ x > i ( T2 - > - > }o, which is invariant under a r- 
indcpcndcnt "local gauge transformation" 4 } x,i( T ) ~ * 4>x,i( T ) + Ac,i> i s wcll-dcfincd, whereas (<f> x ,i(Ti)4>y,i{ T '2))o ls 
not. 
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From Eq.(2.14), the spatial coupling of V Xt i in the Hamiltonian is given as (neglecting the gauge 
field for simplicity) Hy ~ V£iVx+j,i + C -C, and therefore we replace as V x .i — ► (— )^VJV Xt i in the 
action where \x\ = 1(— 1) for even (odd) site. After this replacement, the first term of f/y becomes 

-\ E [^ + „ + cc.]^ |^y Xii | 2 -(rf-i)^|y Xii | 2 , (2.22) 

where = U x .jU x+i ^V x+hl ~ V x>i . 

Let us define the effective action Sep of the Cooper-pair field V Xi i and the gauge field as, 

Z = j [dVdA]e- ScF " s ^ = j [d^dA]e~ s . (2.23) 

Then up to the second-order terms, 

Scp ~ Jdr\\Y,\D j V x A"+ J -^Y,\ D r V M\ 2 

L x,i^j x.i 

-(JAi +d-l)J2 \V x ,i\ 2 + J dk d 2 [e- 2 V(k) - l]- l A {-k)A {k) 



where D T V Xji = (d T - ieA xfi )V Xji . 

At long distances, the last term of S'cp in Eq.(2.24) becomes the usual relativistic kinetic terms 
of A xfi , 

dk k 2 A (-k, T)A (k, t), (2.25) 



I dT I' 



whereas the effects of the short-range part of V x . y generate terms like |T4,i| 4 - Similar quartic terms 
also appear in the higher-order terms of the cumulant expansion. Furthermore, the second and third 
terms of H in Eq.(2.2) are expressed in terms of V x ^ by the perturbative calculation in powers of 
C, D in the cumulant expansion. 

In the derivation of the GL model action S'cp given in this section, only the lowest-order terms 
of 5p Xi i are considered. It is obvious that higher-order terms like 0(((5/9a;^) 2 )-terms in Eq.(2.8) give 
nonvanishing contributions to the GL model. However we think that the full GL model, which 
contains all relevant terms for studying the SC phase transition, has a similar form to that obtained 
in the above. The expansion in powers of the order-parameter field V Xi i, the invariance under a 
7r/2-rotation in the spatial planes and the invariance under a local gauge transformation, 

U x>i ^ e^'+'U^e-*"' , V^^e-^+'V^e-™*, (2.26) 

which originates from Eq.(2.4), restrict the form of the GL action. Effective field-theory model for 
the granular s-wave SC has been derived in a similar way [6]. Then in order to study phase structure, 
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order of phase transitions, etc, we propose action of the GL model as follows, which is a kind of the 
U-V model proposed previously, 

-d 2 Y (U x ,iU x+iJ V x+jti VZ J + C.C.) -c 3 J2 (V^iVUijV^VZj + C.C.) 

x,i x,i 

where (x) = A Xtfl +A x+ ^-A Xj i, — A x +„ ^ (/i, v — 0, 1, d) and we have introduced a lattice also 
in the r-direction. It should be noticed that the gauge field U x ^ (A x ^) also has the r-component 
U x .o = e tA:c ' besides the spatial ones 2 , whereas the Cooper-pair field V x ^ has only spatial compo- 
nents. For numerical calculation, we fix values of the parameters ci, • • • ,/3 in Eq.(2.27) and study 
phase structure by varying g. For large <?, the electromagnetic interactions between V Xj i is getting 
weak, and the hopping of V x< i is enhanced. Because of the uncertainty relation between the number 
n x ^i and phase ifx.i, An^.j • Aip Xt i > 1, a SC state, in which V x /s are stabilized, is formed. For 
small <?, on the other hand, charge density n Xj j tends to have definite distributions and as a result 
an insulating Wigner crystal or the normal metallic state appears. In other words, the parameter g 
plays a role of 1/h 



3 Numerical results in (2+l)D 

Let us turn to the numerical studies of the present lattice U-V model for the SC phase transition 
in two and three spatial dimensions [7]. In order to investigate phase structure, we measured the 
"internal energy" E and the "specific heat" C, 

E = -(S G l)/V, C=((S gl -(S gl )) 2 )/V, (3.1) 

where V = L 3 is the system size. For investigation on 2D thin film SC materials, we consider 3D 
symmetric lattice as we are studying the QPT at T = 0. In the action Sql in Eq.(2.27), we fix 
the values of c±, ■ ■ ■ , [3 and measured E and C as a function of g. As we explained in the previous 
section, the parameter g controls quantum fluctuations, and therefore an observed phase transition 
is nothing but a QPT. 

We mostly studied the cases with the parameters c\ = c 2 = c 3 = — d 2 = 1, a = (3 = 3, 5, 10 and 
also the London limit of V x ^. The negative value of c? 2 means that the configurations corresponding 

It should be noticed that in the action (2.27) V Xt i couples with both U x fi and U x +i t o- This means that V Xt i has 
electric charge at x and x + i as the original Cooper pair given by Eq.(2.1) does. 
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to the d x 2_ y 2-wa.ve SC are enhanced for large g. In all these cases, we have observed no anomalous 
behavior of E, whereas we found that C exhibits a sharp peak as g is varied. We show the system- 
size dependence of C in Fig.l for the cases of a = /3 = 5 and the London limit. It is obvious that 
the peak develops as L is getting larger in both cases. Similar behavior of C was observed for the 
case with a = /3 = 3, 10. Then we conclude that a second-order phase transition takes place as g is 
increased. 




Figure 1: System-size dependence of the "specific heat" C for a — (3 — 5 (left) and the London 
limit (right) of the (2 + 1)D UV model. The peak develops as L is getting larger in each case. The 
result indicates that the second-order phase transition occurs at around g — 0.902(0.490) for the 
a — (3 — 5 (London limit) case. 



In order to verify that the observed phase transition is a SC transition, we measured gauge- 
invariant gauge-boson mass Mq (the inverse of the magnetic penetration depth) . Mq is defined as 
the inverse of the correlation length of the operator 0(x) = Yli - =1 2 e ij sin(i*ij(x)), where e%2 — 
—£2i = 1[8]. In order to evaluate Mq accurately, we introduce Fourier transformed operator of O(x) 
in the 1-2 plane, 

O(x ) = J2 0(x)e vlxl+ip2X2 . (3.2) 

XI, X2 

Then we expect the following behavior, 

(O(s +t)O(x )) oc e -Vpi+P 2 2+ M G\ (3.3) 

In the numerical calculation, we put pi = P2 = \ and verified that the correlator (O(xo + t)0(xo)) 
exhibits the exponential decay as in Eq.(3.3). More precisely, we define Mq as Mq = sign(A 2 — 
p 2 )^/\ 2 — p 2 , where A is the inverse correlation length of the correlator of O(xo). In Fig. 2, we show 
the calculations of Mq for a — (3 — 5 and the London limit cases. The negative value of Mq in Fig. 2 
is the finite-size effect [8]. The results show that in both cases Mq develops continuously from the 
phase transition point obtained by the measurement of C . At first-order phase transition points, the 
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Mq exhibits a sharp jump from vanishing to a finite value[3]. Therefore we conclude that the phase 
transition observed by the measurement of C is a second-order SC transition. We have observed 
similar behavior of Mq for the other cases with a = (3 = 3, 10. 




0-88 0.9 0.92 0.94 q.46 0.47 0.48 0.49 0.5 0.5 i 0.52 

g g 

Figure 2: Gauge-boson mass Mq for the cases of a — j3 = 5 (left) and London limit (right) of the 
(2 + 1)D UV model. Mq develops at the critical coupling obtained by the measurement of C. This 
result indicates that the SC phase transition occurs and the order of the transition is of second order. 



4 Numerical results in (3+l)D 

In this section, we shall show results of the numerical study of the (3+l)D cases. As most of the 
materials of the unconventional SC have quasi-2D structure, studies on 2D and 3D systems are useful 
to obtain the physical picture of the real materials. It is known that quasi-2D systems except for 
thin film samples exhibit 3D properties in critical regions. 

As in the (2 + 1)D case, we first show the calculation of C in Fig. 3 for the case a = (3 = 5. There 
exist two peaks at g ~ 0.67 and 0.78. Then we carefully studied the system-size dependence of each 
peak and the results are shown in Fig. 4. From the calculations in Fig. 4, we conclude that the both 
peaks are signal of second-order phase transition as they have the system-size dependence. Similar 
behavior of C was observed in the case with a = (3 = 10 (see Figs. 5 and 6) and the London limit. 

In order to investigate physical meaning of the phase transitions, we measured the gauge-boson 
mass Mq. Results are given in Figs. 7 and 8. Mq starts to develop at the first critical coupling 
g c <~ 0.67. Furthermore at the second critical coupling g' c <~ 0.78, Mq slightly changes its behavior 
upwards. Therefore, the SC phase transition occurs at g c ~ 0.67. On the other hand at g' c ~ 0.78, 
a certain mechanism suppressing fluctuations of the gauge field A Xtfl further starts to work. 

It is quite helpful to study the monopole density of U and V gauge fields for understanding 
physical meaning of the second phase transition. In the ordinary s-wave SC state, vortices appear as 
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Figure 3: Specific heat C for (3 + 1)D U-V model with a — f3 — 5. System size is 8 4 . There exist 
two peaks at g ~ 0.67 and 0.78. 



x 



0.665 0.67 0.675 0.68 0.685 




Figure 4: System-size dependence of two peaks in C for a = j3 = 5. Result indicates that both peaks 
are signal of second-order phase transition. 




Figure 5: Specific heat C for a = p = 10 with system size 8 4 . There are two peaks. 
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Figure 6: System size dependence of two peaks in C for a = /3 = 10. Result indicates that both 
peaks are signal of second-order phase transition. 
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0.7 0.75 0.8 



Figure 7: Gauge-boson mass Mq in (3 + 1)D system for the a = (3 — 5 case. Mq starts to develop 
at the first critical coupling g c ~ 0.67. Furthermore at the second critical coupling g' c ~ 0.78, Mq 
slightly changes its behavior upwards. 
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Figure 8: Gauge-boson mass Mq in (3 + 1)D system for the a = (3 = 10 case. Mq starts to develop 
at the first critical coupling g c ~ 0.68. Furthermore at the second critical coupling g' c ~ 0.84, Mq 
slightly changes its behavior upwards. 



low-energy excitations. In the GL theory of the s-wave SC, the vector potential A em and the order 
parameter are given as follows for vortex configurations [9], 

A cm = ^-grad 6, * = |*|e ie , (4.1) 

where 

I A cra ■ dl = = -^2nn, n = ±l,±2,---. (4.2) 

7 2e 2e 

In the present U-V model, similar vortex configurations exist. However as the F-field is defined on 
links, the F-field has also monopole-like configurations just like the ordinary compact U(l) gauge 
field, i.e., vortices can terminate at monopole and anti-monopole. The phase transition from the SC 
to normal states occurs as a result of the condensation of (infinitely long) vortices, but we can also 
expect that some phase transition, which is related with monopole dynamics, occurs at a certain 
coupling constant. 

In Fig. 9, we show the calculation of the U and V-monopole densities. We employ the definition 
of the monopole density in the lattice gauge theory given in Ref.[10]. Please remark that the F-field 
is a 3D "vector field" and its monopole configuration can be defined straightforwardly. On the other 
hand, the [/-field is a (3 + 1) D "vector field", and we define the U-monopole just focusing on its 3D 
spatial component. 

From Fig. 9, it is obvious that the U-monopole density pu is very low in the whole parameter 
region of g as the [/-field is the noncompact gauge field. It starts to decrease very rapidly at around 
the first phase transition point g c <~ 0.67. At the second phase transition point g' c <~ 0.78, pu 
is substantially vanishing. On the other hand, the V-monopole density pv starts to decrease at 
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g c ~ 0.67, but it has a finite value between the first and second phase transitions. After the second 
transition at g' c <~ 0.78, pv decreases very rapidly. Similar behavior is observed also in the case of 
a = (3 = 10 and the London limit. (See Fig. 10.) This result indicates that a finite density of short V- 
vortices connecting monopole and anti-monopole survive as low-energy excitations in the parameter 
region between the first and second phase transition points. From the above consideration, it is 
expected that the second phase transition is a very specific one to the U-V model and it does not 
exist in the ordinary s-wave SC. 

0.0175 r 1 1 1 ^ 

0.015 - \ 
0.0125 - \ 
0.01 - \ 

§ \ 
0.0075 » 

0.005 V 

0.0025 - 

o t , 1 • ' i - 

0.6 0.65 0.7 0.75 0.8 

g 

Figure 9: U and V-monopole densities, pu and py, in the (3 + 1)D U-V model with a = (3 = 5. 





Figure 10: U and V-monopole densities, pu and pv, in the (3 + 1)D U-V model with a = (3 = 10. 



5 Discussion 

In this paper, we have derived the effective gauge model of the Cooper-pair field sitting on links 
from the microscopic quantum Hamiltonian by the path-integral methods. This lattice GL theory 
for unconventional SC is a lattice gauge model with dual gauge fields. There the electromagnetic 
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field U x ,n is a noncompact U(l) gauge field, whereas the link Higgs field V x ^ is a compact U(l) 
gauge field. Then we studied phase structure of the model mostly focusing on the QPT at T = by 
means of the MC simulations. For the (2 + l)D case, we found that there exists a second-order phase 
transition for the normal to SC phases as the coupling constant g controlling quantum fluctuations is 
increased. The gauge-invariant gauge-boson mass start to develop at the critical coupling obtained 
by the measurement of the specific heat C. On the other hand for the (3 + 1)D case, we found 
two second-order phase transitions. From the results of the measurement of Mq and the monopole 
densities Pu(v)i we identify the first one as the normal-SC phase transition and the second one as 
the transition of V-monopole suppression. 

It is quite instructive to compare the above results to those obtained in our previous papers, in 
which some related models are studied. In Rcf.[2], we investigated the compact U(l) U-V model in 3 
and (3 + 1)D in which U x ^ is the compact U(l) gauge field. In particular in section 6 of Ref.[2], we 
studied the QPT of the model and found that a first-order phase transition from the confinement 
to Higgs (SC) phases occurs. On the other hand in the present paper, we found two second-order 
phase transitions for the noncompact U-V model in similar parameter regions with the above. This 
result can be understood as follows; in the compact U-V model, the U and U-fields are (almost) 
symmetric though there is no time component of U-field. The monopole proliferation transition of 
the U and U-gauge fields takes place simultaneously as actually observed in Ref.[2]. On the other 
hand in the noncompact case, there is the asymmetry between the U and V gauge fields. Therefore 
roughly speaking, the two second-order phase transitions in the noncompact case coincide in the 
compact case and as a result a first-order phase transition appears. Phenomena similar to the above 
are expected to occur rather generally. In fact in the very recent paper [11], we studied multi- flavor 
Higgs U(l) lattice gauge models in various parameter regions and found that when two second-order 
phase transitions coincide at certain points in the parameter region of the model, a single first-order 
transition appears. One may ask if there exist parameters controlling the distance between the 
two second-order transitions in the present model. We expect that the coefficient C3 in the action 
(2.27) controls the distance, e.g., the two peaks get closer for larger value of c 3 , as the topologically 
nontrivial excitations of the V- field are suppressed for large C3. This problem is under study and 
the result will be reported in a future publication. 

Preliminary result of the numerical study on this problem seems to support the above expectation. 
Furthermore, we found that the relative magnitude of two peaks in C changes as the value of C3 is 
varied. For larger C3, the second peak is getting larger. If the two second-order transitions merge and 
as a result a single first-order phase transition appears as c 3 is getting large, interesting behavior of 
the amplitude of U-ficld v = \V x j\ is expected to be observed. From the behavior of py in Figs. 9 and 
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10, it is expected that for g < g' c the potential term in Sql, —f3v 2 +av 4 7 determines the expectation 
value of v as (v) <~ \J~§^ because the phase degrees of freedom of V x j fluctuate strongly. As g is 
getting larger than g' c , the other terms in Sql start to contribute to (v) and as a result (v) is a sharp 
increasing function of g near g <~ g' c . At the first-order phase transition pint, on the other hand, (v) 
is expected to exhibit a sharp discontinuity and have a hysteresis loop as pv does. In any case, the 
above interpretation of the structure of the phase transitions in the U-V models is useful and also 
interesting. 

In Rcf.[3], we studied thermal phase transition of the noncompact U-V model in 3D. There not 
only the electromagnetic gauge field U x ,i but also the Cooper-pair field V Xti are put on all links of 
the cubic lattice, i.e., the direction index i = 1,2 and 3. We found first-order phase transitions from 
the normal to SC phases in the parameter regions c 2 > 0, d 2 = and c 2 = 0, d 2 < 0. As the results 
obtained in the present paper for c 2 > 0, d 2 < show the existence of the second-order QPT in 3D, 
it is possible that the thermal phase transition becomes of second order in the parameter region like 
c 2 > 0, d 2 < 0, i.e., there exist tricritical points. This problem is under study and results will be 
reported in a future publication. 

Let us turn to the discussion on the QPT of the high-T c cuprates given in the previous work. 
In Ref.[12], the QPT of the cuprates was studied by assuming that the critical behavior can be 
described by a 3D quantum XY model whose Hamiltonian is given by, 

Hxy = ^ ^2 nxV x , y n y - ^ ^ J cos(ip x+i - ip x ), (5.1) 

where h x is the number operator conjugate to the phase of the Cooper pair ip x on the site x of the 
cubic lattice, [n x , ip y ] — iS x . y . Detailed investigation on the model (5.1) was given for the short-range 
interaction case, V x , y = VoS x . y , and an effective GL theory in terms of the order parameter ip x ~ e lVx 
was derived. Near the quantum critical point, the superfluid density p s was obtained as follows by 
means of a mean- field approximation neglecting fluctuations of ij) x , 

Ps = l(^)\jo- V /4d). (5.2) 

The critical point is identified as 




and for ^ < Ad > Ad), the system is in the SC (normal or insulating) phase. Though we 
studied the case of the long-range interactions Eq.(2.3) in this paper, we can qualitatively expect 
the relation like Vb oc 1/g, J cx g. Then the above result concerning to the SC QPT is consistent 
with that obtained in this paper. Furthermore as Mq oc \/l/p s in the SC phase, Eq.(5.2) predicts 
Mo cx \/(Vq/Jo)c — (Vq/Jo) oc — 5c, where g c is the critical coupling constant. The numerical 



1G 



results obtained in this paper qualitatively coincide with this behavior of Mq and are consistent 
with the experiments[13], though in our model the Cooper pair is put on lattice links instead of sites 
and behaves like a compact gauge field. 
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A Derivation of the GL theory 



In this appendix, we shall prove Eqs.(2.11) and (2.13). We first prove Eq.(2.11). To this end, let us 
consider the following quantity, 

H A = /dfe^Ie-^W-ll-^oC-^ioW + ^^^.-eA^^j-e^.o) 

= T^JdkY, {[e~ 2 V(k) - I]" 1 + l}A Q (-k)A (k) - -^i{-k)A {k) + ^M-k)^(k) 

i,j L 

= hl^x* h^ io( - fc) -^^ ( - fc) )( io(fc) -^rW^ (fc) ) 



+v- 1 (k)M-k)^(k) 



where T(fc) = [1 - e^" 1 ^)]" 1 . Therefore 



e- s * = J[dA ]e-f dTHA . 
Next we shall prove Eq.(2.13). To this end, let us consider the following integral, 

= \ E [ J ~ leiA {Ki Je-^e-*-™) (V x+J , Je- iA e^) 



.J e i (<l>x,i-<t>x+j,i-A)_^ cc ^ 



where A = A x+ ij + A x j. Therefore 

e -fdrH t 



J [dV*dV]( 



- f cItH v 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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